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ABSTRACT 

Context. Westerlund 2 is a young, massive, obscured stellar cluster of our Galaxy. It harbors the most massive star with well deter- 
mined parameters, WR20a (82+83M Q ), a dozen very early O-type stars (03-7), and a wealth of PMS stars. Although of importance, 
this cluster is still not well known. 

Aims. The high-energy properties of this cluster, especially those of its early-type stars are examined in details. The variability of the 
X-ray sources is investigated. 

Methods. A monitoring of the field was performed using three Chandra observations. This dataset probes daily as well as monthly to 
yearly timescales and provides the deepest X-ray view of the cluster to date. 

Results. The two Wolf-Rayet stars WR20a (WN6ha+WN6ha) and WR20b (WN6ha) were analyzed in detail. They are both very 
luminous and display very hard spectra, but WR20b does not seem to vary. On the contrary, WR20a, a known eclipsing, colliding- 
wind binary, brightens in the X-ray domain during the eclipses, i.e. when the collision is seen face-on. This can be explained by the 
properties of the wind-wind collision zone, whose high density leads to a large absorbing column (2 x 10 24 cm -2 ). 
All twelve O-type stars previously classified spectroscopically, two eclipsing binaries previously identified and nine newly identified 
O-type star candidates are detected in the high energy domain; ten of them could be analyzed spectroscopically. Four are overlumi- 
nous, but the others present typical L X /L B0L ratios, suggesting that several O-type objects are actually binaries. Variability at the ~ 2cr 
level was detected for a majority of the sources, of unknown origin for the putatively single objects. 

Faint, soft, diffuse emission pervades the entire field-of-view but no clear structure can be identified, even at the position of a blister 
proposed to be at the origin of the TeV source HESS J1023-575. 

Finally, the X-ray properties of PMS objects were also investigated, in particular the brightest flaring ones. They provided an additional 
argument in favor of a large distance (~8kpc) for the cluster. 
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' 1. Introduction of the O-type star population and on the eclipsing characteris- 

^ tics of WR20a. iDame I d2007l) recently discussed the distance of 

. £j 1.1. The Westerlund 2 cluster ^ Westerlund 2 c i ust er on the basis of the CO emission of the 

K*N ' The young (l-2Myr) stellar cluster Westerlund 2 lies in the Hn § iant molecular cloud MC8 (likely associated with the cluster) 

b ! region RCW49. It is an intense star forming region with numer- and the H i 21 cm absorptions along the line of sight towards 

- - ' ous Pre Main Sequence (PMS) object s, as enlighted by recent Westerlund 2. This author found that the cluster must be located 

Spitzer data (IChurchwell et al. L 120051 and references therein). on the far side of the Carina s P iral ™ m and inferred a kinematic 

The cluster also contains many early-t ype stars: at least tw o distance of 6 0+ 10 kpc in coarse agreement with the value of 

Wolf-Rayet stars, WR20a and WR20b (Ivan der HuchtT |200T1) . |Rauw et al. | (|2007|). An indirect argument can be added in favor 

and twelv e very hot O- stars with spectral types ranging from 03 of a ' hi g h ' value, using the properties of the diffuse emission 

to 06.5 dRauw et al. L 12007b . In the IR and radio domains, these detected by |Townsley et al. | <|2004j). The latter authors found an 

stars appear surrounded by two shells, one enclosing WR20a and X " ra y luminosity of about 8-9x10 erg s (for an adopted 

the core of Westerlund 2, the other being around WR20b. distance of 2.3 kpc). Considering that the diffuse emission from 

There are still many uncertainties about the properties and massive star forming regions typically amOTints_to2-5 X 10 32 erg 

content of this cluster. For example, its distance is still a mat- s P er early-type star (j Townsley et al. I |2003i 12004)), we would 

ter of intense debate . To cite only the m ost recent papers: ex P ect values of 3 " 7 x 10 er 8 s for a cluster containing at 

lAscenso et al. I (120071) and iTsuiimoto et all (120071) proposed a least 14 ob J ects earlier than 07 (including the two Wolf-Rayet 

'low' value of 2-5kpc, based on the X-ray and IR luminosi- systems). This is compatible with the observed emission only if 

ties of PMS stars; while iRauw et al. I J2007I) rather argue in fa- the distance is rather large (4-7kpc rather than 2.3kpc). In this 

vor of 8.0 ± 1 .4 kpc, based on the spectrophotometric properties P a P er < we Wl11 thus 8k P c as ^ distance to Westerlund 2. 
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Fig. 1. Three-color images of Westerlund2 in X-rays (a, left: full field; b, right: central 8' in the middle and c, below: core of the 
cluster). The red, green and blue colors correspond to soft (0.3-1. keV), medium (1.-2. keV) and hard (2.-10. keV) photon energies, 
respectively. The images were corrected by the exposure maps and adaptively smoothed prior to RGB combination. 



1.2. X-ray properties 

Studying the X-ray emission of young clusters provides cru- 
cial information, especially if the cluster is heavily obscured 
as is the case of Westerlund2. PMS objects, mainly T Tauri 
stars, are more easily detected in the X-ray and IR domains, 
and their high-energy flares can be analysed to get their decay 
times, which are related to the physical size of the flaring region. 
In the high-mass range, the interacting colliding wind binaries 
and young magnetic objects are more easily identified and their 
properties constrained in the X-ray domain. Finally, the large- 
scale interaction with the surrounding medium is better studied 
through the analysis of the diffuse X-ray emission. 

The X-ray properties of Westerlund 2 have been studied in 
increasing detail with years. The first detection of the cluster 
was made by the Einstein observatory: iGoldwurm et aTl (fl987) 
reported the presence of one X-ray emitter at the position of 
RCW49, interpreted as unresolved sources or truly diffuse emis- 
sion (fro m a supernova remnant or a w ind-blown bubble). Using 
ROS AT. iBeiloni & Mereghettil d 19941) detected possible diffuse 
emission and 24 point sources (7 of them beingin our ACIS- 
I field-of-view), among which HD90273, MSP183 and WR20a. 
The high X-ray luminosity of WR20a was suggested as a possi- 
ble evide nce of an on-go i ng col lidin g wind phenomenon. More 
recently. iTownslev et al. I d2004l) and lTsuiimoto et a!71 d2007l) re- 
ported the results of a first Chandra observation. Thanks to its 
high spatial resolution, this facility is the first to resolve the indi- 
vidual sources in the core of Westerlund 2. More t han 400 X-ray 
sources were identified by Tsuiim oto et al. I d2007l) . mostly PMS 
stars and a few early-type objects. Faint diffuse emission is also 
present; it extends to the south-east of t he cluster and presents a 
soft spectrum plus a very faint hard tail dTownslev et al. 1 12004). 

To investigate further the high-energy properties of 
Westerlund 2, especially those of its massive star population, we 
have obtained additional X-ray observations, which combine to 
provide the deepest X-ray view of the cluster. It enables us to 
study the variability of the X-ray sources on several timescales 
(short-term variations within one exposure; monthly to yearly 
changes between the datasets). 

The paper is organized as follows. Section 2 presents the 
datasets and their reduction. Sections 3 and 4 focus on the O- 
stars and WR-stars content, respectively. The search for a coun- 
terpart to the recently discovered TeV source and the presence of 
diffuse emission are discussed in section 5, while section 6 con- 
centrates on the flaring sources. Finally, we conclude in section 
7. 

2. Observations and Data Reduction 

2.1. Journal of observations 

Westerlund 2 was observed three times with the ACIS-I cam- 
era onboard the Chandra X-ray Observatory, first on 2003/08/23 
for 38ks (PI Garmire), then on 2006/09/05 and 2006/09/28 both 
for 50ks (PI Rauw), see details in Table Q] The first dataset was 
taken in very faint mode, whereas the last two were taken in 
faint mode; in all cases, the frame time was the 3.2s standard. 

1 The MSP numb ers come from the numbering scheme devised by 
Moffa t et al. I jl99lh . 



Table 1. Jo urnal of observatio ns. Phases <p refer to the ephemeris 
of WR20a dRauw et al. Ll2007h . 



Name Start time End time (p 

ObsT 2003/08/23, 18:21:23 2003/08/24, 04:53:56 0.34-0.46 

Obs 2 2006/09/05, 15:52:01 2006/09/06, 06:29:22 0.28-0.44 

Obs 3 2006/09/28,05:46:20 2006/09/28,20:14:43 0.41-0.57 



For all three observations, the ACIS-I chips covered a similar 
field-of-view but the position angle was slightly different from 
one exposure to another. 

The datasets were downloaded from the CXC websiufl and 
we have followed their recommendations for the reduction and 
analysis of the data. As they had been processed by the latest 
version of the CIAO pipeline, we did not completely reprocess 
them. 

2.2. Detailed data processing 

Background flares were only detected for the first observation, 
and times with background count rates larger than 0.02 cts s _1 
were therefore discarded. The astrometry was improved by com- 
paring the position of the bright X-ray point sources WR20a and 
MSP 18 with their counterparts in the NOMAD catalog^. Shifts 
in pixels of (-0.19,0.20), (-0.10,0.04), and (0.14,0.20) were ap- 
plied using the CIAO task wcs_update to the first, second and 
third observation respectively. The event lists were then repro- 
jected without randomization, first to apply this astrometric cor- 
rection, then to match the X,Y coordinates from the last obser- 
vation. A merged event list was created using the CIAO task 
dmmerge (with [subspace -expno] to avoid an uncorrect merg- 
ing of the good time intervals). Source detection was performed 
on this merged event list using the CIAO wavelet algorithm. It 
was applied separately to 1024x1024 images of the whole field 
(pixel of 1.4"), the central 8' (pixel of 0.5") and the core of the 
cluster (pixel of 0.25"). The derived source positions were used 
to define source/background regions for use in estimating count 
rates or in extracting spectra. 

Exposure maps were then calculated for each observation 
and each energy band, assuming a monochromatic source of 
energy 1 keV (total band, 0.3-10. keV), 0.7 keV (soft band S, 
0.3-1. keV), 1.5 keV (medium band M, l.-2.keV) and 3 keV 
(hard band H, 2.-10. keV), respectively. Merged exposure maps 
were obtained using the CIAO task reprojectjmage. The result- 
ing fluxed images (in units photons crrT 2 s px~ 2 ) were used 
to derive background and exposure-corrected photon fluxefl 
The merged data were also used to make the three-color images 
shown in Fig.Q] 

Spectrum extraction was done for each source presenting 
more than 100 counts in the merged event list. It was performed 

2 http://cxc.harvard.edu/ 

3 http://cdsarc.u-strasbg.fr/viz-bin/Cat7I/297 

4 For a typical source with absorbed thermal emission (JV^L West2 = 
10 22 cirT 2 and kT=0.3 keV), a photon flux of 1 photon cirT 2 s corre- 
sponds at the time of our observations to a count rate of 350 cts s" 1 , 
an absorbed flux of 1.7x10"' erg cirT 2 s _I and an unabsorbed flux of 
5.7xl0~ 8 erg cirT 2 s _I . 
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separately on each observation. As recommended by the CXC, 
we extracted the spectral response files with the tasks mkacisrmf 
and mkarf (the latter including a mask file when the considered 
source was close to a gap and affected by the spacecraft's dither- 
ing). When there were less than 100 counts per observation, the 
individual spectra were merged using the FTOOLS tasks math- 
pha, addrmf, addarf. Lightcurve were extracted for the brightest 
sources using the CIAO task dmextract. 

Throughout our analysis, we have used CIAO 3.4, CALDB 
3.3.0.1 andXspec 11.2. 

3. O-type stars 

3. 1 . Detection in the optical range 

The first spectroscopic observations of the early -type popula- 
tion i n Westerlund2 revealed six 06-7 V: stars (iMoffat et al. L 
I1991I) . The presence of such 'late' objects seemed incompatible 
with the presence of a very massive binary such as WR20a (see 
below). A new observing campaign was therefor e undertaken , 
which led to a serious revision of the spectral types dRauw et al. L 
I2007I) : at least twelve O-type objects are present in Westerlund 2, 
with several of them displaying a very early type (0 3^4). All 
twelve O-type stars catalogued by iRauw et aTl (120071) are de- 
tected in the X-ray domain (see Fig. |2), but their characteristics 
differ significantly from one to another (s e e belo w). 

From the photometry of IRauw et al. I (|2007), we could also 
select O-star candidates: they have colors typical of early-type 
stars belonging to the cluster (B-V between 1.3 and 1.6 mag, 
see Fig. 7 of that paper) and absol ute magnitudes small er than 
-3.9, a value typical of 09.5V stars dMartins et al. LI2005I) . Some 
of these candidates are also definite X-ray emitters. The opti- 
cal properties of this subgroup, composed of 9 objects, are pre- 
sented in Table [2] Their color excess was derived from the ob- 
serve d color using an i ntrinsi c color of -0.27, typical of O-type 
stars (Ma rtins & Plez L 12006). Then an estimation of their spec- 
tral type was made from their absolute V magnitude at 8kpc: they 
are late O -type stars, Q7-Q9.5. T he relevant bolometric correc- 
tion from [Martins & Plez I d2006l) was then applied to find an 
estimate of their bolometric luminosities (Table . 



Fig. 2. The core of the Westerlund 2 cluster seen in visible (a, 
left) and X-rays (b, right). X-ray contours superimposed on 
the visible image (c, below; the X-ray data were adaptively 
smoothed). 



Table 2. Properties of candidate late-O stars. The last two lines 
repro duces the data of two eclipsing binaries (from lRauw et al. L 
120071) . 



Name 


RA,DEC (J2000) 


V+ov 


B-V±cr s _i/ 


MSP32 


10:24:03.78,-57:44:39.8 


15.355+0.010 


1.285+0.014 


MSP 120 


10:23:58.44,-57:45:13.0 


16.21 1±0.009 


1.539+0.016 


MSP 165 


10:23:55.18,-57:45:26.9 


15.604±0.009 


1.602+0.013 


MSP168 


10:24:01.61,-57:45:27.8 


14.829±0.031 


1.260+0.039 


MSP441 


10:24:01.45,-57:45:31.3 


14.918±0.091 


1.283+0.138 


Src 1 


10:24:21.28,-57:47:27.6 


14.671±0.022 


1.412+0.025 


Src2 


10:24:16.25,-57:43:43.8 


14.914±0.015 


1.597+0.017 


Src 3 


10:24:06.64,-57:47:15.8 


16.473±0.011 


1.644+0.019 


Src 4 


10:24:02.28,-57:45:35.3 


13.302±0.039 


1.371+0.059 


MSP44 


10:24:00.5,-57:44:45 


15.67 


1.25 


MSP223 


10:23:59.2,-57:45:40 


15.75 


1.50 



binaries detected in the X-ray domain (see Table[3]0. The extrac- 
tion region used for each star was a circle whose radius r, given 
in the second column of the table in Chandra pixels, is a com- 
promise chosen to include most source photons while avoiding 
contamination from nearby sources; background regions were 
surrounding annuli or nearby circles devoid of sources. Errors 
on the photon rates and luminosities were estimated from the 
relative error on the net counts. Hardness ratios were calculated 
from the photon fluxes in the S, M, and H bands using HRj =(M- 
S)/(M+S) and HR 2 =(H-M)/(H+M). The X-ray luminosities in 
the 0.5-10.keV band were evaluated for a distance of 8 kpc us- 
ing spectroscopic fits (only rough ones for the faintest objects) 
and were dereddened only for the interstellar column. 

The brightest objects had enough counts to provide a de- 
tailed X-ray spectrum and they were thus also analyzed spec- 
troscopically. The derived best-fit parameters are shown in 
Table [4] The second column provides the origin of the spec- 
trum: Obs 1, 2, 3 or Mgd (in case the individual spectra did 
not present enough counts and the individual data from the 
3 observations had to be merged prior to analysis). The fit- 
ted model has the form wabs(A^" t ) *wabs(Af H )*mekal(kr) 
or wabs fJV? t ) * [wabs (JVf ) *mekal (kT] ) +wabs (N 2 H ) * mekal(kr 2 )], 
with wabs(A^) fixed to the interstellar column dRauw et al. L 
120071; iBohlinetall 1 19781) . In the table, quoted fluxes and lu- 
minosities are in the 0.5- 10.0 keV energy range and the unab- 
sorbed luminosities L™ abs are corrected only for the interstellar 
absorbing column. Note that lightcurves and spectra were ex- 
tracted in the same regions as defined above. An Lx/Lbol figure 
(Fig. H]l summarizes our findings. 



We have added the 07 star HD90273 to the list of O-type 
objects to analyze. Note that, since it is much more luminous 
and less absorbed than the stars of Westerlund 2, it is probably a 
foreground object, however we have calculated its luminosity as 
if it were an 07 supergiant at 8kpc. 

Spectral types, bolometric luminosities and ab- 
sorbing columns (deri ved from the colo r excess using 
Aff t =5.8xl0 21 xE(B- V). iBohlin et al. I Il978l) for the known 
O-type stars (from IRauw et al. L 12007b and the candidates 
objects (calculated from the photometry as explained above) are 
shown in Table [3] 

3.2. X-ray properties 

We first derived the photon rates using the fluxed images for 
the known O-type stars, the candidates late-O and the eclipsing 



3.2.1 . Bright X-ray emitters 

Since the first detection of massive stars in the X-ray domain, it 
is known that their X-ray luminosity scales with their bolometric 
luminosity followi ng Lx ~ 10~ 7 Lbol (the so-called 'canonical' 
relation, see e.g. 

BerghoferetaO Il997l) . Recent results have 
shown that this relation is rather tight for homogeneous groups 
of stars (NGC6231 s ee ISana et al. 1120061 and Carina OBI see 
Antokhin et al. 2007). Outliers are generally peculiar objects, 
i.e. magnetic stars or colliding-wind binaries. In Westerlund 2, 
four of the bright sources clearly lie above the 'canonical' re- 
lation with Lx/Lbol ratios larger than -6.4. The brightest one, 
MSP18, presents a log(Lx/LgoL) of -5.9 (to be compared with 



5 Note that the values of the photon rates given in that table are raw, 
i.e. uncorrected for the encircled energy fraction, whereas the X-ray 
luminosity, derived from spectral fits, include that correction. 



Table 3. Properties of the early-type stars in the optical and X-ray domains. Background- and exposure-corrected source fluxes (in 10~ 5 photons cirT 2 s -1 ) are estimated for 
each observation from the fluxed images in a region of radius r (in px) and for the 0.3-10.keV band. The last columns give the average hardness ratios and dereddened X-ray 
luminosity in 10 32 erg s and for the 0.5-10.keV energy band (see text for details), a colon (:) indicates an uncertain value. 



Star 


Sp. Type 


log(L SOL 


/As,) C/10 22 r 


Mean ph. flux 


Obs 1 


Obs 2 


Obs 3 


HR, 


HR 2 


/ untibs 
X 


log(Ax7Asoz.) 


WR20a 


WN6ha+WN6ha 


6.28 


1.13 


4 


11.87+0.17 


10.3±0.4 


11.2+0.3 


13.2+0.3 


0.64+0.03 


0.385+0.012 


51.7" 


-6.15 


WR20b 


WN6ha 


5.86 


1.04 


9 


3.77±0.09 


3.69+0.18 


3.72+0.15 


3.88+0.16 


0.61+0.08 


0.656+0.016 


19.0" 


-6.16 


HD90273 


071: 


5.80 


0.22 


20 


0.81±0.06 


1.03+0.11 


0.77+0.09 


0.66+0.10 


-0.83+0.02 


-0.79+0.24 


2.82" 


-6.93 


MSP18 


05.5V-III 


5.89 


0.88 


3 


6.92±0.14 


7.4+0.4 


6.4+0.2 


7.4+0.2 


0.38+0.05 


0.19+0.02 


35.8° 


-5.92 


MSP32^ 


09.5V: 


4.61 


0.90 


3 


0.025±0.014 


0.07+0.04 


0.006+0.022 


0.03+0.02 


-0.15+0.50 


-0.92+1.38 


0.15 


-7.02 


MSP44[ 


09.5V: 


4.44 


0.88 


4 


0.23±0.03 


0.27+0.08 


0.23+0.04 


0.23+0.05 


1. 


0.11+0.12 


0.53 


-6.31 


MSP 120* 


09.5V: 


4.58 


1.00 


4 


0.050±0.018 


0.06+0.05 


0.08+0.03 


0.01+0.02 


-0.07+0.47 


-0.47+0.59 


0.41 


-6.54 


MSP151 


07III 


5.37 


0.95 


2 


0.14±0.02 


0.12+0.04 


0.17+0.03 


0.13+0.03 


0.16+0.29 


-0.28+0.16 


0.57 


-7.19 


MSP157 


05.5V 


5.37 


0.92 


1.5 


0.25±0.02 


0.35+0.05 


0.21+0.04 


0.20+0.04 


0.27+0.20 


-0.10+0.10 


1.16° 


-6.89 


MSP 165* 


07V: 


5.09: 1.08 


4 


0.52±0.04 


0.48+0.06 


0.63+0.09 


0.45+0.08 


0.57+0.16 


-0.49+0.09 


1.98° 


-6.37 


MSP 167 


06III 


5.14 


0.98 


2 


2.05±0.07 


1.80+0.12 


2.25+0.11 


2.00+0.13 


0.57+0.07 


0.00+0.04 


8.85" 


-5.77 


MSP 168* 


08V: 


4.90: 0.89 


2 


0.079±0.017 


0.05+0.03 


0.08+0.02 


0.11+0.02 


-0.23+0.27 


-0.42+0.25 


0.41 


-6.87 


MSP171 


05V 


5.58 


1.00 


4 


0.16±0.03 


0.13+0.03 


0.18+0.04 


0.13+0.09 


0.22+0.27 


-0.93+0.35 


1.09 


-7.13 


MSP175 


05.5V-III 


5.49 


0.91 


2 


0.41±0.03 


0.34+0.05 


0.41+0.05 


0.45+0.05 


0.34+0.14 


-0.56+0.08 


1.78° 


-6.82 


MSP182 


04V-III 


5.45 


0.93 


4 


0.18±0.03 


0.20+0.04 


0.12+0.04 


0.22+0.06 


-0.24+0.16 


-0.76+0.20 


0.58 


-7.27 


MSP183 


04V 


6.07 


1.03 


2 


0.70±0.04 


0.67+0.07 


0.60+0.06 


0.86+0.09 


0.09+0.10 


-0.44+0.06 


4.49" 


-7.00 


MSP188 


04V-III 


5.93 


0.96 


2.5 


3.43±0.08 


3.65+0.17 


3.41+0.14 


3.28+0.14 


0.23+0.05 


-0.10+0.03 


14.4" 


-6.35 


MSP 199 


03V 


5.61 


0.97 


2 


0.11±0.02 


0.10+0.03 


0.10+0.03 


0.12+0.03 


0.50+0.42 


-0.56+0.18 


0.57 


-7.44 


MSP203 


06V-III 


5.84 


0.95 


2 


0.57±0.04 


0.53+0.07 


0.49+0.05 


0.69+0.06 


0.12+0.11 


-0.35+0.07 


3.17" 


-6.92 


MSP223* 


08.5V: 


4.80: 1.00: 2 


0.074±0.015 


0.04+0.03 


0.08+0.02 


0.08+0.02 


0.19+0.41 


0.18+0.21 


0.25 


-6.98 


MSP263 


06.5V 


5.45 


1.09 


4 


0.20±0.02 


0.19+0.04 


0.24+0.04 


0.18+0.03 


0.30+0.18 


-0.47+0.13 


1.17" 


-6.96 


MSP441* 


08V: 


4.89 


0.90 


1.7 


0.084±0.015 


0.09+0.03 


0.07+0.02 


0.10+0.03 


0.25+0.40 


-0.07+0.19 


0.30 


-7.00 


Src 1* 


06V: 


5.31 


0.98 


5 


0.071±0.015 


0.07+0.03 


0.10+0.03 


0.05+0.02 


0.15+0.32 


-0.91+0.34 


0.43 


-7.26 


Src 2* 


08.5m: 


5.22 


1.08 


6 


0.46±0.03 


0.54+0.07 


0.46+0.05 


0.40+0.05 


0.63+0.13 


-0.45+0.08 


1.95" 


-6.52 


Src 3* 


09.5V: 


4.61 


1.10 


4 


0.036±0.010 


0.05+0.02 


0.038+0.017 


0.023+0.015 


> 0.95 


-0.40+0.32 


0.13 


-7.08 


Src 4* 


09.5I: 


5.46 


0.95 


1.5 


0.146+0.018 


0.06+0.03 


0.25+0.04 


0.10+0.03 


1.02+0.21 


0.26+0.12 


0.57 


-7.28 



" fromTableg] 
* OB candidate 
c eclipsing binary 
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Fig. 2. The core of the Westerlund 2 cluster seen in visible (a, left) and X-rays (b, right). X-ray contours superimposed on the visible 
image (c, below; the X-ray data were adaptively smoothed). 

Table 4. Best-fitting models and X-ray fluxes at Earth for the brightest early-type stars in Westerlund 2. For each parameter, the 
lower and upper limits of the 90% confidence interval (derived from the error command under XSPEC) are noted as indices and 
exponents, respectively. The normalisation factors are defined as 10~ 14 \ n e nndV[AnD 2 , where D, n e and «h are respectively the 
distance to the source (8kpc), the electron and proton density of the emitting plasma. Abundances were fixed to solar except for the 
WR stars (for more details, see text). 



Star 


Obs 


N H , 

int 

10 22 cm" 2 


10 22 cm- 2 


kT 
keV 


norm 
10- 3 cm- 5 


Xt (dof) 


fdba 

Jx 

10- 13 ergcirr 2 s~' 


j unabs 

% 

10 33 ergs" 1 


WR20a 
WR20a 


Obs 1 
Obs 2 


1.13 
1.13 


1 53 189 
5 74 14 - 6 

'^1.89 

1 63 1 92 

i.uj 140 


0.87 ? 4 
1 96 6 38 
42° 50 


1 57 16J 

1 -"0.98 

1 07 10 - 7 
4 96 lo:9 

^■ 7U 3.04 


1.15 (33) 


4.51 


4.95 


WR20a 


Obs 3 


1.13 


5 48 7 27 
1 93 224 

3 - 28 240 


1 30 1 52 

35° 43 

1 55 175 


2 61 387 
9 60 2L6 


1.26 (97) 
1.02(119) 


4.10 
5.04 


4.88 
5.69 


WR20b 


Obs 1 


1.04 


1 23 188 


3 62 8 05 


0.26°; 38 


1.24 (24) 


2.06 


1.88 


WR20b 


Obs 2 


1.04 


T e286 

2.18 167 


2 47 313 


41 058 


1.30 (35) 


1.97 


1.77 


WR20b 


Obs 3 


1.04 


1 03 148 


5 21 6 72 


22 28 

"•^0.18 


1.97 (38) 


2.29 


2.05 


MSP18 


Obs 1 


0.88 


qO.18 


3 24 418 

J -^2.53 


q 3 ,0.36 
"• J 0.28 


0.92 (18) 


2.88 


3.70 


MSP18 


Obs 2 


0.88 


qO.13 

u -o. 


3.30 3 * 4 


27 030 

"•■"0.26 


1.10(59) 


2.54 


3.24 


MSP18 
MSP 167 


Obs 3 
Obs 1 


0.88 
0.98 


q0.11 

0.35° 88 


3-04|| 
1.58 107 


32° 36 
1 1° 19 

"• xl 0.08 


1.00 (56) 
1.02 (12) 


2.88 
0.52 


3.78 
0.77 


MSP 167 


Obs 2 


0.98 


Q0.22 


2 32 2 72 


10° 12 


1.95 (21) 


0.72 


1.10 


MSP 167 


Obs 3 


0.98 


0.68 15 


1.28 2 ; 32 


16° 28 

"• lu 0.09 


0.95 (11) 


0.56 


0.78 


MSP188 


Obs 1 


0.96 


45 069 


1 991.40 

1 ' z ' z 1.00 


24° 33 
"•^0.20 


1.56 (26) 


0.88 


1.51 


MSP188 
MSP188 
HD90273 


Obs 2 
Obs 3 
Mgd 


0.96 
0.96 
0.22 


036° 52 
70 84 
0.01JJ 38 


1 33 186 

1 20 

97 1:09 

n 4,0.52 

72 


22 26 
30° 38 
0.02" 17 


0.73 (33) 
1.52 (30) 
0.88 (33) 


0.88 
0.74 
0.16 


1.49 
1.33 
0.28 


MSP157 


Mgd 


0.92 


o °- sl 


1 53 192 


0.010JJ 018 


0.97 (8) 


0.057 


0.12 


MSP 165 


Mgd 


1.08 


1 77 248 


0.28°; 38 


1.66 3 - 2 


0.59 (10) 


0.088 


0.20 


MSP175 


Mgd 


0.91 


1 31 188 
'•■"us 


0.21j$ 


2.25f 8 


1.20(15) 


0.057 


0.18 


MSP183 


Mgd 


1.03 


56° 85 


63° 82 


O.IOJJ" 


0.94 (20) 


0.14 


0.45 


MSP203 


Mgd 


0.95 


0.19JJ 47 


99 1 14 


03 05 
"•"-'0,02 


1.09 (21) 


0.11 


0.32 


MSP263 


Mgd 


1.09 


1 09 1 65 


31° 53 


0.25^ 98 


0.89 (7) 


0.034 


0.12 


Src 2 


Mgd 


1.08 


1 42 171 


39° 68 

"••"0.33 


0.40° 77 


1.20(16) 


0.088 


0.19 



the -6.9 of ISana et al. I l2006h and a high temperature (kr > 
IkeV). The short-term lightcurves of MSP18, however, do not 
present large, significant changes (Fig. [3) but the luminosity 
dropped by 15% in the second observation (Tabl43). Due to the 
limited signal-to-noise ratio, it is rather difficult to judge which 
parameter has changed from one observation to another (flux, 
temperature, absorbing column?), but the high luminosity, the 
hardness of the spectrum, and the variability suggest the star to 
be a colliding-wind binary or a magnetic object. Note that bi- 
narity is also suggested by the position of t he star in the optica l 
colour-magnitude diagram of Westerlund 2 dRauw et al. 1 12007). 
MSP167 and 188 present similar properties, with logiLx/LsoL) 
of -5.8 and -6.3, respectively, and a variability at the 2cr level 
between pointings. Although no large variation are seen for the 
full band short-term lightcurves (Fig. [3), MSP188 exhibits some 
peculiar 2<x change (an increase followed by a decrease) in its 
soft X-ray flux during the last observation. Finally, MSP 165 also 
displays a rather large \og(Lx I Lbol) of -6.4 and some variabil- 
ity (at the 2cr level, the source is 30% more luminous in the sec- 
ond observation), but the source is overall quite soft. All of these 
sources should be considered as potential binaries, and deserve 



a spectroscopic monitoring in the visible domain which is cur- 
rently under way. 

The other bright objects have more typical logiLx/LsoL), 
with values between -6.5 and -7.3. Among these, MSP157 and 
203 display unusually high temperatures. Some of these objects 
also exhibit 2-3cr variations: HD90273 is 36% less luminous in 
the last observation, MSP183 and 203 are 35% more luminous 
in the last observation, the flux of MSP 157 nearly doubles in the 
first observation, whereas that of Src2 increases by 35%. If these 
stars are really single, as their modest luminosity suggests, the 
cause for their variability is hard to define. 

3.2.2. Faint X-ray sources 

The faintest X-ray emitters of this sample (13 out of 26) did not 
have enough counts for a detailed spectroscopic analysis (less 
than 100 counts in total). We must therefore rely only on their 
photon fluxes and rough 'colour' pro perties. In thi s categ ory, we 
find four O-type stars classified by iRauw et al. I d2007l) : from 
their bolometric luminosity and extinction, we would have ex- 
pected at least 300 counts for each of these but we detect less 
than 100 counts. This is especially surprising for MSP199, an 
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Fig. 3. Full band background-corrected lightcurves in units of 
counts per 5ks bin for the three Chandra observations (earliest 
at the top, latest at the bottom). Data from MSP18 are shown in 
the left panel, those of MSP188 in the middle one and those of 
WR20b in the right one. Note that the first and last bins, when 
affected by uncompleteness, were discarded and that the first ob- 
servation of MSP 18 was affected by the spacecraft's dithering 
(thereby explaining the lower number of recorded counts). 

03-4 star that should be at least as luminous as the neighbour- 
ing MSP183 (also 03-4) and certainly more luminous than the 
neighbouring MSP203, an 06 object. However, Fig. [2] clearly 
shows that it is not the case. Only a few of these stars are found 
significantly variable at the 2cr level: MSP 182 is 40% less lu- 
minous in the second observation, MSP32 is nearly undetected 
and Src4 is 3-4 times more luminous in that same exposure; 
MSP168 is 50% less luminous in the first observation, while this 
happens for Srcl in the last one. Overall, their Lx/Lbol follows 
the 'canonical' relation (Fig. |4j, but some of them display large 
HRs (MSP223, Src 4, MSP44). 

In conclusion, most O-type stars appear to follow the Lx - 
Lbol relation derived for the somewhat later O-type stars of 
NGC6231. The larger scatter for Westerlund2 can largely be 
explained by the fainter X-ray fluxes and the poorer knowl- 
edge of the ir properties, nota bly multiplicity, compared e.g. to 
NGC6231 dSana et al. 1 120061) . The situation appears definitely 
more complica ted than the faint+soft vs. bright+hard dichotomy 
suggested by Tsuiimoto et al. (2007, see also Fig. |4j). 

4. Wolf-Rayet stars 

4.1. WR20a 

In 2002, WR20a was serendipitously found to be a binary with 
F=3.7d and e-Q. Follow-up analyses revealed the system to 
be the most massive binary known with accurately determined 
masse s: it consists of two very similar components of 82-83 M 
each (Ra uw et al. L I2005I and references therein). This excep- 
tional object also harbors a wind-wind c ollision that t rigger s line 
profile variability in the visible domain ( Rauw et al. , 2005]). The 
two components have types WN6ha, but they are only 'Wolf- 
Rayet stars in the disguise', i.e. they are likely core-hydrogen 
burning stars with high mass loss rates, only slightly evolved off 
the main sequence. 

4.1.1. X-ray spectrum 

In Ra uw et aTl d2005l) . the wind composition of WR20a was an- 
alyzed: the amount of Helium was found to be consistent with 
solar, but the nitrogen abundance appears six times solar and the 
carbon abundance about or less than 4% solar. We include these 
chemical peculiarities in the spectral fitting, and freeze the abun- 
dances of other elements at solar values. Two-temperature com- 
ponents were necessary to obtain good fits (Table 0). Note that 
the low-temperature component seems the more variable from 
one observation to the other. 

The X-ray spectrum of WR 20a appears overall rather hard. 
This is partially due to the absorption by the rather large inter- 
stellar column density, but reflects also the intrinsic hardness of 
the spectrum. Indeed, our two-temperature fits of the ACIS spec- 
tra of this source indicate temperatures around 1.3-2.0keV for 
the hotter component. We emphasize that plasma at temperatures 
that high is usually not found in the winds of single early-type 



stars, thus suggesting that a significant part of the X-ray emis- 
sion of WR 20a arises in the wind interaction region. Assuming 
that the emission arises from a radiatively cooling shock-heated 
plasma, we find that pre-shock velocities of about 1000 km s~' 
are required to account for this hot plasma. This is twice the 
value of the on-axis pre-shock velocity of the wind- wind interac- 
tion th at we estimated from a simple j3 velocity law dRauw et al. L 
2005). This probably indicates that such a simple approxima- 
tion of the pre-shock wind velocity law is not appropriate in an 
interacting-wind close binary system such as WR 20a. It might 
also mean that the contribution from the off-axis collision, which 
occurs at higher speed, is not negligible. In what follows we shall 
thus assume that the pre-shock velocity amounts to 1000 km s _1 . 

4.1 .2. X-ray lightcurve 

The three Chandra exposures were taken around the secondary 
eclipse: Obsl samples the descending branclQ, Obs2 corre- 
sponds to phases slightly before and at the beginning of the 
eclipse, while Obs3 covers entirely the eclipse (Fig. |5j. A first 
important conclusion from the Chandra light curve of WR 20a 
is the lack of an X-ray eclipse. This result can only be explained 
by a rather large size of the X-ray emitting region. The absence 
of strong variations could be explained if the X-ray emission of 
WR20a were produced in the individual winds, i.e. without a 
strong contribution from the wind interaction. However, there 
are two arguments against this assumption: first the rather high 
plasma temperature (see above) and second the X-ray luminos- 
ity that corresponds to an Lxlh, \ ratio well in excess of that 
of single early-type (actually O-type) stars. Such properties are 
indeed more typical of colliding wind binaries (either O+O or 
O+WR). In the framework of a model where the X-ray emission 
arises mid-way between the stars as a result of a wind interac- 
tion, the absence of an eclipse sets stringent constraints on the 
minimum size of the X-ray emitting reg ion. For an orbital incli- 
nation of 74.5° (see iRauw et al. I I2007 and references therein), 
we find that the lack of an occultation implies a strict minimum 
radius of the emitting region of 12 R©. 

In fact, the lightcurve not only reveals the absence of eclipse, 
it also indicates a brightening of the emission near conjunction 
phase, i.e. when the collision region is seen face-on. More pre- 
cisely, on the fluxed images, WR20a appears 20-30% more lu- 
minous in the last observation while the photon fluxes in the first 
two exposures are compatible within 2-3cr (Table [3]). During 
each exposure, increasing trends are detected, but individual bins 
are generally within one or two sigma of a constant value. The 
last two observations, during which WR20a was not affected 
by dithering, follow each other in phase, and the behaviour of 
WR20a can thus more easily be quantified when the lightcurves 
are folded in phase (Fig. [5j. At the beginning, WR20a presented 
about 175 cnts per 5ks-bin, while at the end it reached 250 
cnts per 5ks-bin. If we compare the lightcurve in the hard and 
medium band^l (see Fig. |3J, the increase appears larger in the 
former, thereby leading to a slightly smaller hardness ratio in the 
last observation. The binary is thus significantly X-ray brighter 
and somewhat softer during the photometric eclipse. 



6 Note that we need to be rather cautious about the results from the 
first observation, since during this exposure WR20a moved in and out 
CCD2 due to the spacecraft's dithering. 

7 Since WR20a is heavily absorbed, there are not enough counts to 
make a lightcurve in the soft band. 
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Fig. 4. Left: X-ray luminosity vs bolometric luminosity. The solid line corresponds to the canonical relation of San a et al71 ((2006); 
the stars represents the Wolf-Rayet objects, the filled circles the brightest O-type stars (see Table |4]i and the open circles the faintest 
ones (see Table[3j. Right: Hardness ratios, symbols as in the left panel. 



Fig. 5. Left: Visible lightcurve of WR20a (top figure), X-ray lightcurves with bins of 2ks (middle figure) and 5ks (for the full energy 
band, bottom figure). Right: Lightcurve in the medium (M, top) and hard (H, middle) bands, together with evolution of the hardness 
ratio HR2 — (H - M/(H + M) with phase (bottom), all shown in 5ks bins. The X-ray data from the second observation are shown 
in solid black circles and those from the third one in open red triangles; uncomplete bins were discarded. The uncertainty on the 
phase cr^ is about .01, taking into account the uncertainty on the period and the time elapsed since our optical observations (~700d, 
iRauw et all 120071) . 



4.1.3. Interpretation 

A possible explanation for this puzzling result could come from 
the orbital modulation of the absorption in the dense wind inter- 
action itself. Since the winds of the two components of WR 20a 
most likely have the same strength, the 'stagnation' point, where 
the hottest plasma and hence the hardest X-ray emission is ex- 
pected, should be located roughly at the centre of mass of the 
binary. The column density to the centre of mass is determined 
mostly by the density of the postshock material and by the thick- 
ness of the postshock region along the line of sight. Near quadra- 
ture phases, the column density along the line of sight is largest 
because the line of sight crosses the material near the centre of 
ma ss where the density is largest. 

lAntokhin et al. I (I2004I) provided a semi-analytical model for 
a radiative wind interaction zoneff These authors propose to cal- 
culate the cooling length Iq as lo/r =s 0.01 v 5 r sin 3 6/M where 
r is the distance between the centre of the star and the stagna- 
tion point in R , v is the local wind velocity in 10 3 kms -1 , 6 is 
the angle between the tangent to the interaction front and the lo- 
cal wind direction and M is the mass-loss rate in 10~ 6 M o yr _1 . 
In the part of the interaction region around the binary axis, that 
is of most relevance for X-ray emission, the deflection due to 
the orbital motion is relatively small and hence 8 - n/2. For 
WR20a, we further adopt r = a/ 2 = 27 Ra, v = 10 00 kms'" 



and M = 8.5 x lO^Moyr -1 (see IRauw etaTll2005l and here- 
above for v). This then leads to Iq 0.86 R Q . Due to the steep 
temperature dependence of the radiative cooli ng, the width of 
the po stshock cooling layer is of order 21q/7 (lAntokhin et al. L 
120041) . which then yields about 0.25 R Q . Antokhin et al. further 
showed that the surface density of the shock region on the binary 
axis <tq should be equal to <x = M/4nva. For WR20a, this 
number amounts to about 0. 1 1 g cirT 2 , implying a mean density 
in the postshock cooling layer of about 6.7 x 10~ 12 gcirT 3 . If 
we assume that the wind interaction region within a radius of 
12 R around the binary axis is delimited by two parallel planar 
shocks, integrating the density along the direction perpendicular 
to the binary axis (using equation [32] of lAntokhin et al. | [2004) 
and a ssuming the same chemical composition as IRauw et aT71 
(2005) yields a (highly ionized) hydrogen column density N H 
of order 2 x 10 24 cm 2 . This extremely large valu^l over (prob- 
ably) only a limited fraction of the interaction region shows that 
the absorption of the X-ray emission by the material of the in- 
teraction region must indeed be quite substantial around quadra- 
ture phases. Of course, this large value of the column density 
only applies to a single point (the stagnation point on the binary 



8 The^ cooling parameter of Stevens et a]7| dl992h amounts to about 
0.02 for WR20a. 

9 Note that due to the steep dependence on the pre-shock velocity, 
adopting v = 500kms~' would imply a column density about 60 times 
larger, and a pre-shock velocity of 1500 km s -1 a value ten times lower. 



axis). Since the X-ray emission arises over an extended part of 
the wind interaction region, the actual (observable) absorption 
will be lower. In addition, in a highly radiative wind interaction, 
various kinds of instabil ities are expected to affect the shape of 
the interaction region (IStevens et al. L 119921) . making it deviate 
significantly from a simple planar region. The latter effect will 
reduce the peak value of the column density towards the stagna- 
tion point, but will also broaden the episode of high absorption 
around quadrature phases. 

4.2. WR20b 



Little is known on WR20b. van der Hucht 1 (120011) assigned it 
a type WN6ha, whereas ISharaetal.l d 199 lb prefer WN7. As 
types and ages are similar for both Wolf-Rayet stars, we have 
used the same chemical composition for the spectral fits (Table 
|4]i as well as the same intrinsic color and bolometric correction 
as for WR20a. Combined with our photometry (V=13.51±0.02, 
B=14.98±0.02), this leads to a bolometric luminosity of about 
7.3xlO 5 L , i.e. slightly fai nter than the faintest component of 
WR20a dRauw et al. L 12001) . 

The photon rates and hardness ratio of WR20b are presented 
in Table [3] while its spectral properties for each observation can 
be found in Table [4] Overall, the system is 37-38% less lumi- 
nous than WR20a both in the visible and in X-rays. The spectral 
fits reveal a large extinction, typical of Wolf-Rayet objects, but 
also a high temperature. The hardness ratio is also large, even 
larger than for WR20a (Fig. 0J. Such properties are only ex- 
pected for interacting binaries. However, the X-ray flux appears 
rather constant on both short and long timescales (Fig. [3] and 
TableO, suggesting the star to be single. This contradictory re- 
sult is important in the context of the current debate about the 
level of intrinsic X-ray emission from WN stars: while som e 
apparently single WN stars, like WR40 dGosset et al. L 120051) . 
remain undetected in X-rays, other objects, e.g. WR6 and 110 
(ISkinner et al. 1120021) . clearly emit X-rays. WR20b seems to be- 
long to the latter category. An optical monitoring is under way 
to better constrain its nature. 



5. HESS J1 023-575 and the diffuse emission 



Aharonia n et al. I ((2007) recently reported on the discovery of a 
Te V high-energy y-ray emission in the vicinity of Westerlund 2 
with the HESS Cherenkov telescope array. This emission ap- 
pears extended with a Gaussian fit yielding cr — (0.18 + 0.02)°, 
equivalent to a diameter of 28 pc for a distance of 8 kpc. The 
source centre lies 40s (or 5.5') west of WR20a. Therefore, it is 
unlikely to be associated with WR20a alone, especially since 
the photon-photon absorptions in the wind interaction zone of 
the colliding wind binary are likely to be extremely severe, pre- 
venting a putative TeV emission from escaping from the in- 
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ner regions of the wind interaction zone of WR 20a. The size 
of the source is actually more typical of a wind-blown bubble. 
lAharonian et al. I d2007[) accordingly propose the TeV source to 
be associated with a blis ter-like structure detected in the radio 
range west of the cluster (Whiteo ak & Uchida 111 9971) . 

To explain the TeV emission, iBednarekl d2007l) proposed a 
hybrid model where both leptons and hadrons contribute to the 
y-ray emission from very young open clusters. In this model, 
the leptons are responsible for X-ray and GeV y-rays through 
synchrotron and inverse Compton radiation, whilst the hadrons 
produce the TeV y-ray emission through n pion decay. The elec- 
trons and the hadrons would be accelerated at the shocks result- 
ing from the interaction of the stellar winds of the massive stars 
with the surrounding matter. Whilst the relativistic electrons lose 
their energy inside the lower density cavity created by the winds 
of the massive stars, the relativistic protons travel over larger dis- 
tances and will only lose their energy when interacting with the 
dense ambient giant molecular cloud. Hence the TeV emission 
produced by these particles should be more extended than the 
X-ray and GeV emission. 

An alternative expl anation has been proposed by 
lAnchordoqui et alTI d2007l) who suggest that the TeV emission 
seen towards Westerlund2 could be the result of Lorentz- 
boosted MeV y-rays emitted in the de-excitation of daughter 
nuclei produced by the photo-disintegration of PeV cosmic-ray 
nuclei in the radiation field of the early-type stars inside the 
cluster core. 

In some of the scenarios discussed hereabove, the TeV emis- 
sion is expected to have an extended counterpart in the X-ray 
domain. However, in the Chandra observation, the diffuse emis- 
sion pervades most of the field and appears rather featureless: 
our X-ray data provide no clear indication of a blister-like struc- 
ture (Fig. [TJ. The brightest emission is to be found to the south- 
east of the cluster, not towards the radio blister, and no bright 
X-ray point source is detected in the close vicinity of HESS 
J1023-575. Moreover, the X-ray spectrum of the diffuse emis- 
sion appears mainly thermal although a faint, hard compo nent 
might (or might not) be non-thermal (Towns lev et al. L 2004). 



6.1. Flaring region 

The spectra and light curves were extracted for the flare observa- 
tion over circular regions with a radius of at least 1 .25 arcsec. For 
isolated sources, the background was evaluated over a source- 
free annulus around the source region. However, in the crowded 
inner regions of the cluster, we rather evaluated the background 
over a circular, source-free, region as close as possible to the 
source. For the lightcurves, whenever possible, we have deter- 
mined the 1/e decay times of the flares (r b s ) by a least square fit 
of the exponential decay (see Table[5]l. Due to the limited num- 
ber of counts, the light curves had to be binned into 5 ks bins, 
thus limiting our ability to determine decay times shorter than 
this binning time. However, all but one source actually have 
r bs > 5 ksec. On the other hand, the spectra were fitted with 
an absorbed (wabs) single temperature, solar abundance, raekal 
model. Unfortunately, the sources are still too faint to extract 
time-resolved spectra over the duration of the flare, and the spec- 
tral analysis was therefore done for the whole exposure. The 
results of the fits are provided in Table[5] The best fit plasma 
temperatures are all quite high (kT > 2.2 keV) as expected for 
pre-main sequence (PMS) stars during a flare. We note that the 
best fit column densities are larger than the interstellar column 
(~ 1 x 10 22 crrT 2 ) expected from the average reddening towards 
the massive stars of Westerlund 2. 

Correlating the positions o f the flaring sour ces with the 
2MASS point source catalogue dCutri et al. U2003I) . we find that 
sources A, B, D, and G have a near-IR counterpart within less 
than one arcsecond. However, the only one with good qual- 
ity 2MASS photometry is Src A (J = 15.47, H = 14.80 and 
K s = 14.46), the other three have highly uncertain measure- 
ments. 

Assuming that the flaring plasma is confined in a closed 
coronal loop as is the case for so lar flares, we can u se the simple 
analytical approach proposed by S erio et alTI d 1 99 ll for an appli- 
cation of the method see e.g. Fav ata et al. 1120051) to express the 
loop half-length I as a function of the maximum temperature at 
the top of the loop r max and the thermodynamic decay time t±: 



6. Variable X-ray sources 

A photometric analysis was carried out on the 0.3-10.keV fluxed 
images of each observation using the standard procedures for 
aperture photometry available under iraf. The positions of the 
sources have been derived from the merged data in order to in- 
crease the signal-to-noise ratio. Apertures with radii between 1 .5 
to 5"were used. The relative errors were estimated from the pho- 
ton statistics derived from the fluxes in the uncorrected images. 

To identify flaring sources, we have compared the photon 
rates of each source in the three observations. Many objects 
appear variable (about 6% at the 99% significance level), but 
we decided to focus on those that display a pronounced inter- 
observation variability (x 2 > 100) and are bright enough to allow 
for a detailed analysis of their light curve and spectrum. In this 
way, we have identified eight objects that display a strong flare 
during one of the observations (see Fig. [6]). It should be empha- 
sized that we are not aiming at establishing a complete census 
of variable sources, but rather at characterizing the properties of 
the brightest flaring candidates. Note that in some cases we find 
that the observation does not cover the entire flare but only the 
decay phase. 



/ ~ 0.132 Tfc 

Here I is expressed in R Q , kT in keV and t± in ksec. When the 
temperatures are determined from fits of ACIS CCD spectra, a 
correction has to be applied to the observed temper ature to derive 
the m aximum temperature at the top of the loop dFavata et al. L 
2005): T mlLX — 0.068 T^ 2i) where the temperatures are in K. The 
observed decay time of the flare exceeds the intrinsic thermody- 
namical decay time due to the effect of heating that can extend 
into the decay phase of the flare. The relation between the two, 
F(£) = T bs/Tth, can be determined from the slope f of the flare 
decay in a log T - log y/EM diagram dFavata et al. L l2005h . As 
mentioned above, the flaring sources discussed here are too faint 
to allow a time-resolved spectral analysis. Therefore, we can 
only determine lxF{£). The corresponding values of lxF(0 are 
3.0 - 4.8 R Q (see Table|5]i. If we assume that the flaring loops are 
freely decaying with no heating, F(£) ^ 2 dFavata et al. L [2005) 
which actually represents the lower limit on F{£ ). Therefore, the 
typical loop sizes are < 1.5 - 2.4 R . The flaring regions of the 
PMS stars in Westerlund 2 are thus rather compact and typical 
of many PMS as well as active main-sequence stars that display 
loops anchored on the stellar photosphere with no interaction 
with an accretion disk. 
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Fig. 6. Light curves of the eight sources that display a strong flare during one of the Chandra observations of Westerlund2. The 
sources are identified as in Tableland the number of the observation where the flare occured is also indicated. For each panel the 
x-axis shows the time (in ksec) elapsed since the beginning of the observation. The light curves were binned with a time bin of 
5 ksec but note that the first and last bins might be affected by uncompleteness). 



Table 5. Spectral analysis of the flaring PMS candidates, ordered by RA. Columns 2 and 3 yield the coordinates (J2000) of the 
sources, whilst the fourth column indicates the number of the observation during which the flare occured. Columns 5-9 list the best 
fit parameters of the spectral fits with a wabs x mekal model. For each parameter, the lower and upper limits of the 90% confidence 
interval are noted as indices and exponents, respectively. Finally, column 10 yields the 1/e decay time of the flare; kr max is the 
maximum plasma temperature during the flare (derived from the fitted kr) and / is the loop half-length. 



Source 
# 


RA 


DEC 


Obs. 
# 


N 

(10 22 cm" 2 ) 


(keV) 


Jx 

(io- 


j~unabs 
14 X 9 1 

erg cm s ) 


Xt (dof) 


Tobs 

(ksec) 


(keV) 


lxF(0 
(Rg) 


A 


10:23:34.27 


-57:48:09.3 


2 


<0.21 


4 05 993 

2.02 


1.4 


1.4 


1.75 (4) 


2.4 






B 

C 


10:23:57.20 
10:24:00.00 


-57:42:59.5 
-57:44:21.5 


1 
1 


9 QQ4.78 

1 35 211 


5 24 79 - 9 
> 3.18 


8.8 
4.7 


16.1 
6.2 


0.39 (8) 
0.47 (9) 


8.7 
6.6 


12.8 


4.1 


D 
E 
F 


10:24:00.86 
10:24:01.58 
10:24:02.05 


-57:45:56.8 
-57:45:22.7 
-57:44:33.0 


3 
2 
1 


2 07 2 - 84 
1 56 219 

1,JU I.I)7 

i -27O.3O 
l.J/ 574 


S Q-342.9 
J ' 7J 3.02 

q q i9.4l 
J ' y 2.30 

> 1.78 


5.2 
6.0 
2.6 


8.9 
10.9 
3.4 


1.12 (11) 
0.74 (16) 
0.76(1) 


6.2 
12.0: 


14.9 
9.0 


3.2 
4.8 


G 


10:24:03.21 


-57:45:34.5 


2 


1 70 2 - 77 

L - /u 0.88 


t 99 15.0 


3.0 


6.1 


1.66 (8) 


10.6 


7.2 


3.7 


Ft 


10:24:10.15 


-57:46:35.9 


3 


3 00 507 


2 17 9 ' 72 


1.2 


4.2 


0.28 (5) 


10.7 


4.5 


3.0 



6.2. Further evidence for a large distance 

Tsui imo to et al. I ((2007) derived a loose constraint on the dis- 
tance towards the cluster of 2 - 5 kpc, based on the median X-ray 
luminosity of PMS stars in the mass range 2.0 - 2.7 M . In fact, 
it seems rather unlikely that the X-ray data can provide stringent 
constraints on the distance. If we consider for instance the X-ray 
luminosity of the flaring sources, we find that their mean, ab- 
sorption corrected, Lx amounts to 6.7 x 10 32 ergs-' during the 
flares if we adopt a distance of 8.0kpc. While this is a rather high 
X-ray lumino sity, we note that f or the brightest flaring sources 
in NGC 6231, ISana et al. I J2007h derived L Y ~ 3 x 10 32 ergs _1 . 
In addition, the brightest flaring sources in the Orion Nebula 
Cluster observed within th e COUP project disp lay peak lumi- 
nosities of 8 x 10 32 ergs- 1 (Fava ta et al. L 12005). The luminosi- 
ties of our flaring sources in Westerlund 2 are thus not untypical, 
especially since they are probably biased towards higher values: 
first, we have selected only the brightest flaring objects; second, 
these luminosities were corrected for the full column densities 
derived from the spectral fits. As noted above, these column den- 
sities are larger than the values estimated from the interstellar 
reddening of the early-type stars in the cluster core and it could 
thus be that our absorption-corrected X-ray fluxes are somewhat 
overestimated. 

Our X-ray data might however still give hints for constrain- 
ing the distance of Westerlund 2. Indeed, a total of 202 X- 
ray sources have a n optical counterp art at less than 1" in the 
B V photometry of Rau w et aTl (12007) (see Fig. [7j». The colour- 
magnitude diagram of these sources indicates that most of the 
X-ray sources with an optical counterpart appear along the main 
sequence of the cluster or are foreground field stars. The latter 
objects actually provide a large fr action of the counterpar ts, a 
fact apparently underestimated by Tsuiimoto et al. I (120071) . In 
addition, the colour-magnitude diagram shows no sign of the 



Fig. 7. (B—V, V) colour-magnitude diagram of the X-ray sources 
that have an optical counterpart within a correlation radius of 
less than 1". The dashed lines correspond to the main-sequence 
between spectral types 03 and B2 for a distance modulus of 
14.52 and reddened with A v = 4.68 (leftmost curve) or 5.84 
(rightmost curve) and assuming Ry = 3.1. 



ZAMS turn-on. Indeed, none of our flaring source^] seems as- 
sociated with an optical source in our photometry. The optical 
counterparts of these sources must thus be fainter th an V = 21 . 
For v ery young clusters such as e.g. NGC 6231 dSana et al. L 
|2007|) . the optically brightest PMS stars that display X-ray emis- 
sion are found to have absolute magnitudes of My - 0.5 - 1.5, 
in good agreement with the predictions of evolutionary models 
dSiess etaU 120001) . For a distance of 8. kpc (DM = 14.52 ) and 
interstellar extinction A v = 4.7-5.8 (Ra uwet aL~l,l2007l) . the 
optically brightest PMS stars in Westerlund 2 would be expected 
to occur somewhere in the range V = 19.7 - 21.9. Were the clus- 
ter at 2kpc, the turn-on would be brighter by 3 mag and would be 
easily spotted in our data, and counterparts to our sources could 
be detected. Their absence therefore suggests a large distance for 
the cluster. 



7. Conclusion 

In this paper, we have analyzed three Chandra observations 
aimed at monitoring the massive objects of Westerlund 2. In this 
heavily obscured region, 24 O-type stars (identified or candi- 
dates) were detected and analyzed. Most of them display typical 



10 Note that source A falls outside the field-of-view covered by our 
photometric data. 
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values of their Lx/Lbol ratio, though the scatter is rather large 
(about 60%) due to the rather low signal-to-noise ratio of the X- 
ray spectrum, but also the much poorer knowledge of their bolo- 
metric luminosities and multiplicities. The most X-ray luminous 
objects (MSP18, 165, 167 and 188) are probably binaries and 
need a specific follow-up. About half of the O-type stars also ex- 
hibit variability at the 2<x level or more, whose origin is unknown 
in these supposedly single objects. 

WR20a, the most massive binary, displays a hard and heav- 
ily absorbed X-ray spectrum. It brightens by about 40% in the 
X-ray domain during the photometric eclipses. The X-ray emis- 
sion is believed to come mostly from the wind-wind collision 
region, which would easily explain the presence of a high tem- 
perature in the best-fit models. If the colliding wind region is 
sufficiently large and opaque, a modulation will appear when the 
system rotates and it should then be brighter when seen 'face- 
on', as is the case during eclipses. The other Wolf-Rayet star of 
the field, WR20b, is thought to be single and its flux is accord- 
ingly constant. However, the star also exhibits a high tempera- 
ture and a high X-ray luminosity, and its hardness ratio is even 
larger than for WR20a. These latter characteristics are more typ- 
ical of colliding-wind binaries. 

The photometry of the other X-ray sources was examined 
for each observation. Many objects appear variable, and we fo- 
cused on the eight brightest, flaring ones. Their lightcurve decay 
time and spectral properties enabled us to infer rather compact 
flaring regions. The comparison between the optical and X-ray 
properties of the X-ray sources again suggest a large distance for 
Westerlund2. 
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